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Abstract 

The qualitative results for string rotation in the frame of Relativis- 
tic Flux Tube Model and the quantitative and qualitative results for 
decay of massive open string states of string theory are used as the ba- 
sic principles in the Monte Carlo implementation for decay of massive 
open string. The presented model can be used as an ingredient into 
any Monte Carlo model of multiparticle production for different types 
of collising particles. In the interval of total cm. energy from 3 GeV 
{e~^e~ annihilation) to 1800 GeV (proton-antiproton interaction) the 
presented "soft" (i.e. withaut qq,gg,qg scattering ) model has agree- 
ment with experimental data on transverse momentum distribution of 
secondary particles up to 4 GeV. According to widely known theoretical 
hypothesis, the secondary particles from interval pT >lGeV are results 
of hard scattering QCD states. Therefore (maybe) one can recover a 
question about relation between the classical string solutions of string 
theory and solutions of QCD equations. 
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1 Introduction 



If we remove the quarks from QCD, there would remain a nontrivial SU{3) 
Yang-Mills theory which must have its own spectrum of states. These states 
will become the "glueball" states of QCD (See Ref.|jl|]) . The simple quark 
model can be extended to incorporate such gluonic degree of freedom. In the 
bag model, such approach is to proceed by analogy with the "constituent 
quark" to posit the existence of a " constituent gluon" with the quantum num- 
bers (J^*" = 1 ) of a gluon of weak-coupling perturbation theory 0. In the 
non-relativistic Flux Tube model |jl|] by Isgur and Paton there is the long-range 
strong-coupling limit in which the gluonic degree of freedom have condensed 
into collective string like flux tubes. A flux tube is a directed element (or 
"string"), and quark(antiquark) acts as a unit source of the three-dimentional 
flux tube. Therefore a "junction" can annihilate (be created) there , i.e. there 
are flux tube breaking and fusion. Such theoretical picture leads (see ) to 
an understanding of the linear confinement (quark-antiquark) potential, to lin- 
eary rising Regge trajectiries Hence, in the potential flux tube model, 

a simple QCD-motivated potential |Q, ^ V{r) = + a-r (here k ~ 0.5 and 
a ~ O^GeV"^) is used for calculation of interaction along to string (flux tube) 
lenght. 

In the frame of classical analysis Olson, Olsson and Williams have 
demonstrated that for large angular momenta the leading relativestic QCD 
corrections can be interpreted as the angular momentum and angular energy 
of a rotating flux tube. Therefore , the simplest interpretation of these rela- 
tivistic corrections is to consider the momentum as well as the energy of the 
interacting field. So, the relativistic flux tube model 0] , arises as the 
successor to the potential model. 

In the Relativistic Flux Tube (RFT) model (see Fig.l) the energy of 
motion of the i-th quark (fermionic) is Hi = a-p^i + (3-mi. The tube is inserted 
by the four momentum substitution |Q p^' — > — Pt > where p^ is the new 
canonical four momentum and pt is the tube four momentum computed by 
integration along the tube ^ 

{Pti)r = , 

{Pu)e = Lu = ^{Hu - a ■ r.(l - ^L)°-^), (1) 
Hti = a - Tii^u ■ Ti sin^9^i)"\ 
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where dj^i is transverse velocity of i-th. quark. As one can see , the tube 
substitution is analogous to the introduction of the four potential in QED. 
For small quark masses mi -C ar,, it can be shown that d±i —>■ 1 and the 
quarks are dynamically unimportant for large angular momenta. In this limit 
we recover (see ) the Nambu-Goto string with Regge slope = 

The basic assumption of the RFT model ( see ) is that the QCD dynam- 
ical ground state for large quark separation consists of a rigid straight tubelike 
color flux configuration connecting the quarks. The similar physical picture is 
for classical string solutions in Ref. |10| where Mitchell, Sundborg and Turok 
had considered the decay of massive open string states on the leading Regge 
trajectory. These correspond to classical string solutions in the form of rigidly 
rotating rods whose ends move at the speed of light (see above limit 'd±i — ^ 1 ). 
For these states the lenght L of the string is proportional to the mass {L oc M) 
. In [|r^ Mitchell, Sundborg and Turok had found that the decay rate for a 
string of a "classical" lenght L is proportional to L^*^"^"^)/^^ oc M^'^~^^^^^'^ , 
where c? is a space dimension. For example, for our case d=4 and the decay 
rate F oc L-o-^^ ^ m-°-^\ 

Hence, rotation of string playes the fundamental role in both the flux tube 
model (which has correct theoretical basis for low and middle mass calcula- 
tions) and the string theory which has theoretical description for decay up to 
high mass of string (see ).The question about rotation of string has deep 
physical connection with a question about relation between radial and orbital 
excitations, i.e. relation between momenta Pr and Pg in Fig.l . 



As it had described in Refs.[ll] , |12| for the case of large orbital exci- 



tations L ^ rij. (or Lz ^ Ur , where is the radial quantum number) the 
system behaves as the transverse rotating string (i.e. for the quaziclassical 
long distance approximation it is Pg > Pr in Fig.l ). It is an analogy with so 



called dominantly orbital excited states |T^. In the opposite case L -C (i.e. 
Pg < Pr in Fig.l ) the string is nearly pure inert and almost does't contribute 



into the kinetic part of the Hamiltonian [IT . 

For example, in the Lund fragmentation scheme [0 , for the force field 
between quark and antiquark, an approximation of quasiclassical relativistic 
massless string is realizatedQ. Hence this string does not lie on the Regge tra- 
jectory, has not orbital momentum and rotation (i.e. the momentum Pg in 
Fig.l is neglected). As it is emphasized in [|13| the surface spanned by the 

^The string model becomes semiclassical when we account for its breaking [l6|. 
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massless relativestic string is always a minimal surface. This means that (see 
||13||), in the Lund fragmentation scheme , all interior properties of massless 
relativistic string is determined by the boundary. Therefore in this approx- 
imation the string is nearly pure inert and only (anti) quark as the excited 
end-point of string has energy-momentum. Therefore, the momentum Pg is 
neglected and the process is one-dimensionalQ. For high radial excitation of 
end-point, this approximation is correct to a high degree of accuracy. 

At the some time , after Lund string fragmentational procedure, the final 



state consists of a set of string pieces [|13[, which have the energies, momenta 



and masses, i.e. ones lie on the linear Regge trajectories and have rotation at 
rest. Hence, in Fig.l, the value Pg can't be neglected. Therefore in the Lund 
fragmentation model there are the classical massive (rotating) string solutions 
only for small string masses (meson and baryon hadron states). 

The lund approximation is starting point for our model. According to 



results of Ref. |10| in our string fragmentation scheme we suppose existence of 



the classical massive (rotating) string solutions not only for small mass string 
states, but for the mean and large mass string states too. 

In the Relativistic Flux Tube model and for classical string solu- 

tions of string theory |T^, for large quark separation a rigid straight tube like 
color flux configuration connects the quarks. At the some time, in the soft 
(i.e. withaut qq,gg,qg scattering ) proton- (anti)proton interaction, and in the 
proton-electron interaction, and in electron-positron annihilation for high cm. 
energy , the above both quark and antiquark have long distance strong inter- 
action , ultrarelativistic opposite momenta (in their center of mass system) 
and , due to orbital momentum conservation, they move with non-zero im- 
pact parameter. Therefore for the short time period between formation and 
decay of string , there are physical conditions (for formation and decay) which 
are the same as in the flux tube model or in the classical string solutions of 
string theory [0. Hence, at high energy interaction, the cascade decay of 
massive rotating open string can be serious realistic candidate for description 
of multiparticle production by the use of cascade multistring production. 

The primary string (arising in electron-positron annihilation , proton-electron 
or proton- (anti)proton interaction) has radial direction which is directed along 

^In the Relativistic Flux Tube model the process is three-dimensional, because the 
quark acts as a unit source of the three-dimentional flux tube (string) which has the mass, 
momentum-energy and agular momentum and lies on the Regge trajectory. 
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to summary momentum vector of both the quark and antiquark which act 
as units sources of the three-dimensional primary string (flux tube). There- 
fore there are correspondingly two degree of freedom for excitation of meson 
(baryon) system in its rest frame : the values Pe and Pr (in Fig.l). In Refs. 



[]T^,|T^ it is shown that for conventional meson states, where the flux tube is 
in its ground state, the transverse confinement leads to Gauss law distribution 
for transverse degree of freedom and the string transverse size is about a few 
hundred MeV (i.e. ~0.3 GeV). We use this theoretical result in our model, 
where the longitudinal momentum Pr grows with the mass of string. Therefore 
in our model there is the average summary physical picture : 

-the string is breaking along to direction which is determined by angle 6; 

-the angle 6 is determined by relation between the momenta Pg and Pr 
from the previous decays; i.e. rotation of daughter string has direct 
connection with rotation of mother string 

-the momentum Pg has "restricted" value (which is determined by 
Gauss law) , whereas Pr grows with the string mass. Therefore the string 
of very large mass is breaking longitudinally (i.e. the angle 6 is about 
zero), and the string of small mass is breaking approximately isotropically; 

-the average value of angle 6 grows with number of decays of string. This 
is natural result , which corresponds to classical (rotating) string solutions 
and to the Relativistic Flux Tube model where the string rotation 



is one in a few basic principles. 



As it is noted above, in Ref. |T0[ , Mitchell, Sundborg and Turok considered 



the decay of massive open string states and found that the decay rate for 
a string of "classical" length L (and of mass M) is proportional to L~^'^^ , 
i.e. r oc M~°'^^ for the space-time dimension d=A. The branching ratios 
to differ mass level states illustrate this behaviour graphically in Fig.l (in 
Ref.[0), where the dominant decay modes are into one small string and one 
large string. The same physical picture can be arise in our model , if we for 
simplicity suppose^ that the partial decay rate Fj in the differential kinematical 
interval of daughter strings is proportional to the kinematical characteristics 

■^Of course, it is rough additional supposition which (in future MC implementation ) 
should be changed over correct theoretical result. 
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of daughter strings from this interval, i.e. Fj oc {si + Pzi)~^'^^ , where Ei 
and Pzi are energy and ^-projection (i.e. along to Pr momentum in Fig.l) of 
momentum of i-th daughter string in rest frame of decaying string. Therefore 
the (density of) probability for decay of massive open string with mass M 
in to daughter strings with above kinematical characteristics from differential 
interval is = ^ oc ^"^'^j^fo.sa — = ( ^'J^" )~°'^^ = Zi~^'^^ , where convenient 
variable Zi = ^^^^ is used. Both the zi- and Z2-distribution functions and 
Gauss law for projections of transverse momentum (see above) of daughter 
string in rest frame of decaying string allow to us to determine the masses, 
energies and momentum projections of two daughter strings in the rest frame 
of decaying string (see Fig. 1)0. 

The above theoretical principles and results we use in the algorithm (see 
Sect. 2) which is an attemption to construct the Monte Carlo implementation 
for breaking of massive open string according to quantitative results for decay 
of classical string solutions of string theory and to the results of RFT model 
also. 

The results of Monte Carlo calculations is discussed in Sects. 3 and 4. 



2 Decay of string 

Monte Carlo (MC) algorithm of decay of string is described in this Section. 
The diagram of cascade breaking of string is shown in Fig. 2. The double 
line marks the diquark. The cascade breaking of string is described by the 
iteration of an elementary process presented in Fig. 3 : the string is breaking 
into two non-interacting colourless daughter strings. Three schemes of decay of 
the string into two strings (Fig. 3) are taking into consideration in our model. 
According to the Flux Tube model (see Fig.l) and in the frame of cascade 
logic we generate the break of each string in its rest frame. By every break 
, for example of secondary string of mass M (see Fig. 2) , we know following 
its characteristics (in the Kq c.m.s.0) : the scheme of the decay of secondary 

^According to Relativistic Flux Tube model , and classical string solutions of string 
theory [|o|, for large quark separation a rigid straight tube like color flux configuration 
connects the quarks. Therefore it is necessary to stress, at Monte Carlo simulation of decay 
of above rigid straight tube (string), the zi and Z2 can't be generated independently. 

^It can be, for example, c.m.s. of pp-collisions, or e+e^ annihilation cm. system with z 
axis along to jet axis. 
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string, the flavours of tlie partons participating in tliis decay, tlie mass M, tlie 
energy Eq and tfie momentum vector Pq of the decaying string. 

We generate the break of each string in its rest frame. For transition 
from Kq to string rest frame, let us define new Kl frame, where longitudinal 
momentum of decaying string is equal to zero. This Ki frame moves with 
j3zo = Pzo/Eq velocity along to Zq axis in the Kq frame and with Lorentz-factor 
7o = Eo/M^_ = Eo/{El - P^X^ ' where = (M^ + P]_X^ is transverse 
mass of i-th string. 

In the Kl frame, the decaying i — th string has transverse momentum 
Pio with projections P^o and Pyo. We turn the Kl frame (in the transverse 
xy plane) to obtain the xl axis along to the vector P^q. The parameters of 
this transformation are cosip± = Pxo/Pi.o , smip± = Pyo/P±o ■ Let us define 
this new system as K± one , and let the K± frame move with j3±o = P±o/M± 
velocity and with Lorentz-factor 7_|_o = E±q/M = M±/M , where Ej_q is energy 
of decaying string in Kl frame. Thereby we obtain the K' rest frame of the 
string. 

The longitudinal z' axis of K' frame is parallel to the longitudinal zq axis 
of Kq frame. But according to the Flux Tube model (see Fig.l) the breaking 
string has rotation , i.e. there is nonzero angle 6 between the flux tube and the 
z' axis in our K' rest frame (it is CM. frame in Fig.l). At present we can't 
calculate 6 angle (according to theory and simultaneously ) in frame of our MC 
model ( it is interesting problem which should be overcomed for future MC 
implementation ) , and so for continuity of rotation (of rigid straight tube) from 
the mother to daughter string we determine the polar angle 6 by the formula 
sin6 = P_Lo/Pi(2) ; where Pi{2) is the modulus of the momentum vector of 
first (second) daughter string]^ in the K' rest frame[] , and the azimuthal angle 
if is generated uniformly^ in the interval (0, 2tt). The above formula takes into 
account the 6 fluctuations at fixed P±o and the average 6 angle grows with Pj_q 
of string. In the string rest frame K' the angles 6, ip are determine the line of 
string stretch. Along this line the string break is generated, and so we have 
to determine the new K rest system which is given by the 9, Lp angles in K' 

^The P\{i) value is given by eqs.(3),(20), and so if Flo > -^1(2) i we can use an approxi- 
mation sin Q = \. 

^As it is emphasized in Sect.l , in the string breaking models jl^ , |^ the angle Q (in 
Fig.l) equals to zero. 

^In future MC implementation it is necessary to take into account that the azimuthal 
isotropy can be broken by angular momentum conservation. 
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system and with the y axis along to the vector product^] z ^ z' . 

A simple model for string breaking involves quark pair creation by the 
strong chromoelectric field inside the string (see [jl6l) . The fiavours of quarks, 



produced from vacuum at the decay of string, are generated according to the 
relation 

M : : s = 3 : 3 : 1 . (2) 



In the RFT model IT^IjIIII] the transverse confinement leads to Gauss law 
distribution for transverse degree of freedom and the string transverse size is 
about a few hundred MeV. We use this theoretical result in our model, where 
in the K frame P^i and Pyi components of the transverse momentum of 
first daughter string are generated according to Gauss law (see Sect.l) 

ri{r2) = (y'\2Ti)-^-'' r^'^'''"\xp{-x^/2a^) dx , (3) 

J —oo 

where is uniformly distributed random number in the range (0,1) and a is 
parameter, which is aqq = 0.35 GeV for the string decays in Figs. 3b and 3c. 
For the string decay in Figs. 3a the model parameter a is cru(rf)=0.25 GeV (if 
valence quarks are u(d)-quarks), and the model parameter a is (Js=0.30 GeV 
(if only one valence quark is s-quark). 

In the K frame we determine (see Sect.l) the following variables: 

zt = {Ei + P,i)/M, z+ = {E2 + P.2)/M, z^ = {Ei-P,i)/M, 

z^ = {E2 - P.2)/M, (4) 

where Ei and E2 {Pzi and Pz2) are energies ( the momentum projections 
(on z axis)) of first and second daughter strings. Generally speaking, the 
^-distribution can be singular one . In this case the z-variable should be 
generated from interval ( Zmin, Zmax )• The possible choice for Zmin and Zmax is 
the following. According to law of the conservation of energy, at the production 
of first daughter string we generate the variable zf in the interval 

^Imin < Z]^ < ^imaa; , (5) 

where 

4_ = {E, - \PTr\)lM , (6) 



'Therefore third nonphysical Eulerian angle -0 is absent in our calculations. 
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4ma. = (El + \Prr\)/M , (7) 

where 

El = (M^ + (M{"™)2 - (M2™")2)/(2M) , (8) 

I p;^'^^ I = ({Elf - (M™")2)0-5 , (9) 

where M{"*" and M^*" are the minimal masses of first and second daughter 
strings (see Fig.3). Any daughter string is converted into a hadron if it cannot 
be broken into two hghtest hadrons , i.e. the minimal mass of first (second) 
daughter string equals to the sum of the masses of two lightest hadrons, which 
have the flavours depending on the quark contents of first (second) daughter 
string. Therefore (for simplification) the minimal mass of first (second) daugh- 
ter string can be determined as the double mass of lightest hadron with quark 
composition of first (second) daughter string (see Appendix 1). 

For production of first daughter string, the distribution density of vari- 
able (for all schemes of the string breaking in Fig.3) is parameterized by the 
form (see Sect.l) 

f{z) = z-'-'' . (10) 

The same scaling|^ function (10) can be used for production of second daughter 
string if we generate the variable (see eq.(3)). In the K frame the equalities 
(4) lead to ones : 

zt + z^ = 1 , (11) 

zt . z^ = mi? + {P±f)/M' , (12) 

4 . z^ = {{M,r + {P^f)/M' , (13) 

where Mi and M2 are the masses of first and second daughter strings (see 
Fig. 2). Therefore (after generation of the variable z^ ) the variable 2;^ is gen- 
erated in the K frame (according to the scaling function (10)) in the interval 

^2min < ^2 ^ ^2max ' (-'-4) 

where 

^2-™n = ((Mr-)' + {Pa.?) I {M\l - zt)) , (15) 



^"According to theory of decay of massive open string pO| , in our Monte Carlo implemen- 
tation the distribution of z does not depend on the mass of mother (decaying) string, i.e. 
this distribution is scaling one. 
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Z2ma. = 1 - mr'^Y+{P^f)/{M'zt) (16) 

are determinated according to law of the conservation of energy. In the K 
frame from the conservation of energy and momentum after generation of Zi 
and Z2 we can find (see (4)) the variables 

Zi = I - Z2 , (17) 

4 = 1- zt . (18) 

The energies and longitudinal momenta of daughter strings are given by for- 
mulas (see Appendix 2) 

= M{zt + z-)l1 , (19) 

= M{zt - zi)l2 , 2 = 1,2. (20) 

Because of an equality Zi = {1 — 4) and from eq.(lO) one can see that 
the variable z^ has z'*^'^^ distribution for first daughter string and (according 
to energy conservation) (1 — z)~'^'^^ distribution for second daughter string. 
Therefore (in the frame of presented model) for the fitted values of the free 
parameters, for small mass of mother string and for the random numeration 
of daughter strings, at the decay of mother string into two daughter strings 
the distributions of first and second daughter strings on the part of energy(and 
momentum) of mother string are similar to momentum distributions for valence 
quark and antiquark in mesorQ i.e. similar to f{z) = z~^'^{l — z)~^'^ (see 
Ref. [H ). 

In the models [|T9[,|2y] there is no such result, because for the force field 



between quark and antiquark, an approximation of quasiclassical relativistic 



massless string is realizated (see [|13| and Sect.l) , therefore the scaling func- 
tions are used for the other process (see Fig. 4a and Sect.l). For example, in 
Ref.[Q the scaling function for conversion of the quark into meson+quark is 
equal to /(2) = l-a+3a(l-2)2 , with a=0.88 , in Ref.ll f{z) = (l+c)(l-z)^ 
, with c ^ 0.3 ^ 0.5. 

After determination of the momenta of daughter strings, we control the 
inequalities 

Ei > Pi , 1 = 1,2 (21) 



For z variable there is the same result because the and variables have the same 
parametrization (10). 
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If inequalities (21) are fulfilled for the i-th string, we determine the string 
mass Mi = {{EiY — {PifY'^ and the decay scheme of i-th string and flavours 
of quarks, produced in this decay from vacuum, are generated. Then we control 
the inequality 

Mi > M™^" + M^'"" , (22) 

where the M^^^ is minimum mass of the j-th hadron with the fixed quark 
composition, which can be produced by decay of i-th string. The decay scheme 



(see Fig. 3) is generated according to the relation |]T9 



a : h = 0.95 : 0.05 . (23) 

A model has six parameters ( see eqs. (2), (3), (10), (23) ), which have correct 
theoretical validity. 

If the inequalities (21), (22) are fulfilled for both daughter strings, their 
decays are simulated. The algorithm for decay of daughter string of the mass 
Mi{M2) is similar to algorithm for decay of secondary string of the mass M . 

If the inequalities (21), (22) are fulfilled for only one daughter string (for 
definiteness we take the first string), the decay of the mother string into first 
daughter string and hadron is generated. If < P2 we select the lightest 
hadron from hadrons with the given quark composition. The early generated 
zf , Pxi and Pyi values are not changed, but we generate the hadron mass 
M2 instead of Z2 , if the hadron is resonance. After determination of Z2 from 
(18), we determine Z2 according to (see (13)) 

Z2 = + {Pi)') I (M'zt) , (24) 

and z^ from (17). The energies and longitudinal momenta of the daughter 
string and hadron are calculated from (19), (20) (see Appendix 2). If the in- 
equalities (21), (22) are fulfilled for the daughter string, the algorithm of it 
decay is similar to algorithm for decay of secondary string of the mass M . 

The decay of mother string into two hadrons is generated, if at the decay of 
the mother string into two strings the inequalities (21), (22) are not fulfilled for 
the both strings, or if at the decay of mother string into a string and a hadron 
the inequalities (21), (22) are not fulfilled for the daughter string. There are 
some decay modes for the given quark content of these hadrons. We attribute 
a weight to each decay mode. This weight is equal to the product of three 
factors. The spin factor is equal to (2Ji + 1)(2J2 + 1), where Ji and J2 are 
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spins of the hadrons. The kinematic factor is equal to the two body phase 
space volume or to zero, if the string mass M is smaller than the sum of the 
masses of daughter hadrons. The S't/3-factor is taken into consideration, if 
there are several hadrons with the same quark content, spin and parity. For 
example, Sf/a-factor of 77-meson which is formed from uw-pair is equal to 1/6 
, and the same Sf/s-factor of ti?- meson is equal to 1/2. 

The Ml and M2 masses of resonances are generated after the generation 
of the decay mode. For example, if string decays into two resonances and the 
mass of first resonance must be generated at first, the Mi and M2 masses can 
be generated according to the Breit-Wigner distributions in intervals 

M™" < Ml < M-M^'"" , (25) 

M2™" < M2 < M-Mi , (26) 

where M™*" is maximum sum of masses of particles produced by the decay of i- 
th resonance. In the K frame the square of transverse momentum of resonance 
= P^i + Pyi is generated in the interval (0, P^), and the azimuthal angle is 
generated uniformly in the interval (0, 27r). The momentum of first resonance is 
supposed to have the sharp angle with z axis. The decays of unstable hadrons 
into stable and quasistable particles are generated in the K frame, for example, 
according to algorithm ||2^ . 



The momenta of stable and quasistable particles are transformed from the 
K frame to the K' frame according to formulas 

p'xi = —Pxi cosO cosip — Pyi sirup — p^i sinO cosip , (27) 

P'zi = Pxi sine - p^i cosO , (28) 

v'yi = ~Pxi cosOsirvp + Pyi cosip — p^i sinO simp , (29) 

where p^i , Pyi , Pzi are the momentum projections of z-th particle (in the K 
frame), and 9,ip angles are determined above. After Lorentz transformation 
of the energies and momentum vectors of particles from the K' frame to the 
K± frame, momentum vectors of particles are transformed from the K± to 
Kl frame according to the parameters of transformation cos ip±_ , sin Lp±_ (see 
above) in the transverse xy plane. After Lorentz transformation of the energies 
and momentum vectors of stable and quasistable particles from the frame 
to the Kq frame we know all characteristics of secondary particles in the Kq 
c.m.s. 
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3 e+e annihilation 

3.1 e^e~ annihilation at low energy 



In Sect. 2 a cascade model of string breaking (see Figs. 2, 3) had been con- 
structed as a new model for hadronization. In e^e~ annihilation with center 
of mass energy Ecm. up to a few GeV the process of gluon emission |2^ - 



pO| practically is absent - |3^ . Therefore we deal with only one primary 
quark-antiquark string (see Fig.5). Therefore in this energy interval there 
is unique possibility to check MC cascade model of string breaking (Sect. 2). 
However, our model incorporates the production and decay the hadrons which 
are constructed by u, d, s quarks only. Therefore there is a possibility to check 
this model at center of mass energy only which is under threshold of charm 



quark pair production, i.e. at c.m.energy equals to 3 GeV [31|, |32 



3.2 The formation of the primary string and its decay 

The and e~ annihilate to form a virtual photon which produces a quark- 
antiquark pair (Fig.5). In the presented model the flavour of the quark q 
(and antiquark q ) (Fig.5) is generated according to relation (between the 
probabilities for uu , dd , ss quark-antiquark pair production) 

u: d: s = A: 1 : 1 . (30) 

In the process e+e" qq (in the e~^e~ center mass system) the angle distri- 
bution for the quark is given by the form 

N-^dN/dn = {1 + cos'^e) 3/16TT . (31) 

where 6 is polar angle between the quark momentum vector and electron mo- 
mentum vector, dQ = 2TTsin6d6 . 

A quark q and antiquark q (Fig.5) stretch the primary string A which 
decays into secondary hadrons according to the algorithm of cascade of string 
breaking which is described in detailes in Sect. 2. In e~^e~ c.m.s. the quark 
momentum vector determinates the direction of the string breaking (see Sects. 
1,2). 



3.3 Comparison with experimental data for e+e annihilation 
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In e~^e~ center of mass system (see Fig. 5) the quark momentum vector 
is opposite to momentum vector of antiquark. Therefore after hadronization 
there are two jets of secondary hadrons . A direction for jet axis (see Fig. 6) is 
determinated by Monte Carlo method ( ||3^ ) . 

At Ec.m.= 3 GeV the time for Monte Carlo simulation (see above) of one 
e~^e~ annihilation is 0.40 second. Illustration of a hadronic event from e~^e~ 
annihilation is in Fig. 6. The results of calculations according to the model are 
the curves in Fig. 7 (the multiplicity distribution), Fig. 8 (rapidity distribution), 
Figs. 9 (p_L-distribution). The variables x and x\\ (see Figs.8,9) are defined by 
X = 2P/Ec.m. and = 2P||/Ee,m. (see Fig.6 and Refs. Q, |3 ). 



4 Conclusion 

The last few years have witnessed rapid progress in the description of string- 
string interactions pSf and semiclassical decays of strings (flux tubes). For 
example, in Ref . |16| Gupta and Rosenzweig explored the implications of string 
breaking (flux tube fission) for hadron decays ; and in Ref. Nussinov shown 
that the flux tube model can be successfully applied to analytical calculations 
for various stages of chromoelectric flux tube intersection, rearrangement, and 
eventual fragmentation ("hadronization") of the resulting highly excited flux 
tubes ; in the Ref. [|10[] Mitchell et. all calculated the decay rates for arbitrarily 



massive states on the leading Regge trajectory for open strings. 

In our Monte Carlo (MC) model |^ we attempt to construct the string 
(and its decay) which has the physical properties the same as the string in the 
Relativistic Flux Tube model 0,0)0 and simultaneously the same as the 
classical massive open string solutions [|l^ of string theory. Therefore in our 



MC model the string decay is computed taking the effect of string rotation 
into account, it is the flrst. The second, in the Ref. 0] Mitchell, Sundborg 



and Turok have calculated the decay of the fundamental string in an arbitrary 
number of dimensions. For calculation of decay probability we use decay rates 



I0| for the space-time dimension d=A. 

We can make a comparison with some results on fundamental (Nambu- 
Goto) strings and RFT model. In the decay of a fundamental string |jlO| in 
the critical dimention d=26 Mitchell, Sundborg and Turok flnd that the string 
is likely to break symmetrically (in contrast to the noncritical dimention) into 
two states, each with roughly half the (mass) excitation number of the original 
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state. The same result is for semiclassical decay of excited string states in the 
flux tube model |jl6| , where the space-time dimension (1=4. At the same time 



at the decay of a fundamental string |10| in noncritical dimension the string 



decays by emitting a state close to the ground state and a highly excited state. 
In our Monte Carlo model we use decay rates (from Ref. [|1^]) for the space- 
time dimension (i=4. Hovewer, because of momentum-energy conservation 
we generate the kinematical characteristics of the both daughter strings from 
intervals which are determined by the conservation laws, the first. The second, 
because of isotropy of time-space, and because of P (space inversion) invariance 
of strong interaction, we can suppose that in the statistical set of decays of 
mother string (with fixed mass) the distributions for kinematical characteristics 
for first daughter string are the same as ones for second daughter strin^^ inside 
of the above kinematical interval^. Therefore our MC calculations yield result 
that the share of string decays into the state close to the hadron state and a 
highly excited state is about 25% only. 

The some other properties of our model one can see by comparison with 
electron-positron experimental data at the total center of mass energy equal 
to 3 GeV , where the process of gluon emission practically is absent and we 
can check model at center of mass energy which is under threshold of charm 
quark pair production^. 

The most interesting quantitative results for proton-(anti)proton interac- 
tions will be given elsewher^^, but now it is necessary to stress one interesting 
property which has direct connection with massive open string dynamics |TD 



Because of the rotation (Fig.l) of decaying string (Fig. 3) , the average 
transverse momentum of secondary strings (Fig. 2) grows with the mass of pri- 
mary string ( the A string in Figs. 2, 5). A quantitatively illustration of this 
is in Fig. 10. In other side the average mass of primary string grows with 
the total energy of coUising particles ( e~^e~ , pp , pp ) . Therefore the av- 
erage transverse momentum of secondary strings grows with the total energy 
of collising particles, and so for the interval of total energy (in the center of 
mass system of collising particles) from 3 GeV (e^e~ annihilation) to 1800 

^^At the random numeration of both strmgs. 

^^See choice for intervals and distributions for simulation of zi and Z2 variables from 
Sect. 2. 

^■^In our hadronization scheme only lightest u- , d- , s-quarks is used. 
^^The preliminary results see, for example, in [El|. 
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GeV (proton- (anti)proton interactions[^ there is agreement between experi- 
mental data and our theoretical MC calculations on transverse momentum p± 
distribution of secondary particles up to pt = ^GeV (see Figs.11-13) . It is 
interesting because there is widely known theoretical hypothesis that in the 
p_L interval (approximately) from IGeV to AGeV the particle production is a 
result of QCD hard {qq, gg, qg) scattering states (see |2^). But in our model, 
all secondary particles are result only of decay |]lOl of states of massive open 



strings (classical string solutions of string theory |10| ) and with the same be- 
haviour as for particles which are result of QCD hard scattering statesj^ (at 
least in the p± scale). Therefore (maybe) one can recover a question about 
relation between the classical string solutions of string theory and solutions of 
QCD equations. 

Our results are a hint that it would be very interesting to construct the 
Monte Carlo implementation for the decay of massive open string most closely 
to results of string theoryQ. We hope that our model is first step to this aim. 
But already the presented string breaking model can be used as an ingredient 
into any Monte Carlo model of multiparticle production for different types of 
collising particles ( pp , e^e~ , e~p , heavy-ion reactions ets). 

For heavy-ion reactions there is widely known theoretical prediction that 
hard parton-parton interaction and formation of quark-gluon plasma are two 
sources of secondary particles with high transverse momenta. We hope that our 
calculations shown that the string rotation is also the source of high transverse 
momenta, and so this result can't be ignored in the realistic calculations of high 
energy reactions. 



"'^^In proton- (anti)proton interactions for formation of primary strings, the well-known soft 
version of dual parton model [p2| is used, i.e. hard interactions (gg, gg, qg) between partons 
of collising hadrons is not calculated. 

^^The detailes for production of secondary particles which are the result of QCD hard [qg, 
ffff,^ gg^ scattering states , and the detail discussion about distribution are in Refs. , 

For example, for MC implementation it is best to have the strict (without our additional 
suppositions) theoretical partial decay rate , which depends on the kinematical characteris- 
tics of daughter string, it is the first. The second, for MC scheme for long-distance interaction 
between quark and antiquark, it it necessary to have a theoretical (share or probability for) 
relation between the states with angular (orbital) momentum L — Q (here the quark degrees 
of freedom dominate) and the states with angular (orbital) momentum i 3> 1 (here the 
string degrees of freedom dominate). 
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Appendix 1 



According to experimental data ||3^ for high energy hadron-hadron colli- 



sions, for the produced primary hadron from every hadron multiplet the share 
of the members with big mass grows with total energy of collision. Therefore 
there is second possibility for determination of the minimal mass (of first (sec- 
ond) daughter string), i.e. in second case it is the product of two factors. The 
first factor equals to 2. The second factor grows with energy. For example, 
at \fs < 100 GeV the second factor is the mass of lightest hadron with quark 
composition of first (second) daughter string; at 100 GeV < y/s < 1000 GeV 
the second factor is (with equal probability) the mass of any hadron with quark 
composition of first (second) daughter string; at y/s > 1000 GeV the second 
factor is the bigest mass of the hadron with quark composition of first (sec- 
ond) daughter string. This determination of the minimal mass of the daughter 
strings leads to the better agreement with the experimental data, but it is not 
a decisive factor. 



Appendix 2 

At decay of string of the big mass the sharp angle 6 (in Fig.l) leads to 
the sharp angle between the z axis and the momentum vector of the daughter 
string. Therefore it is necessary to rewrite the formulas (19) as the following 

P,i = M\ zt - z^ I /2 , P,2 = -M I 4 - Z2 I /2 . 



17 



References 

[1] N. Isgur and J. Paton, Phys. Rev. D31 , 2910 (1985). 

[2] J. Donoghue, K. Johnson,and B. Li, Phys. Lett. 99B, 416(1981). 

[3] L. Jaffe and K. Johnson , Phys.Lett 60B, 201 (1976). 

[4] M. G. Olsson , MAD/PH/793 , September 1993. 

[5] E. Eichten, K. Gottfrid, K. D. Lane, and T. M. Yan, Phys.Rev.D 17, 3090 
(1987); 21, 203 (1980). 

[6] Cohin Olson, M. G. Olsson and Ken Wilhams Phys.Rev.D 45, 4307 
(1992). 

[7] M. G. Olsson and Ken Wilhams , MAD/PH/736 , December 1992. 

M. G. Olsson and Ken Wilhams, Phys.Rev.D 48, 417 (1993). 

[9] M. G. Olsson and Sinisa Veseli, MAD/PH/816 , January 1994. 

[10] David Mitcheh, Bo Sundborg and Neil Turok, Nucl.Phys. B335, 621 
(1990) 

[11] A. Yu. Dubin, A. B. Kaidalov and Yu. A. Simonov, |hep-ph/93rT51^ , 
|hep-ph/ 94082121 



[12] D. La Course and M. G. Olsson , Phes.Rev. D 39, 2751 (1989). 
[13] B. Andersson and A. Nilsson , LU TP 92-25 , October 1992. 



[14] M. G. Olsson, |hep-ph/9611340 . 

[15] V.Innocente et al., Phys.Lett., 169B , 285 (1986) 

[16] K. S. Gupta , C. Rosenzweig, Phys.Rev.D 50 , 3368 (1994). 

[17] Frank E. Close and Philip R. Page , |hep-ph/941130Tl 

[18] Cheuk-Yin Wong, Ren-Chuan Wang and Jian-Shi Wu, piep-ph/9412357 



18 



[19] B.Andersson, G.Gustafson, G.Ingelman, T.Sjostrand, Phys.Rep. 97 , 31 

(1983) 

[20] P. Aurenche, F.W. Bopp, J. Ranft , Z.Phys.C 23 , 67 (1984) 

[21] V.Lugovoi , PTI 173-97-HEP, March 1997; 174-97-HEP, June 1997 . 

[22] A. Capella, U. Sukhatme, Chung I Tan, J. Tran Thanh Van, Phys.Lett. 
81B, 68 (1979); A. Capella, U. Sukhatme, J. Tran Thanh Van, Z.Phys.C 
3, 329 (1980) ; A. Capella, J. Tran Thanh Van, Z.Phys. CIO, 249 (1981); 
Phys.Lett. 114B, 450 (1982); A. Capella, U. Sukhatme, C-I Tan, J. Tran 
Thanh Van, Orsay preprint LPTHE 92-38. 

[23] K. Werner, |nucl-th/ 94 11033 . 

[24] B. Andersson, G. Gustafson, Hong Pi, Z.Phys.C 57, 485 (1993) . 

[25] K.Hahn and J.Ranft, Phys.Rev.D 41 , 1463 (1990) . 

[26] V.M.Chudakov and V. V.Lugovoi , Z.Phys.C 59, 511 (1993) . 

[27] T.D.Gottshalk, Nucl.Phys. B214, 201 (1983) ; ibid. B227, 413 (1983) ; 
T.D.Gottshalk,D.A.Morris, ibid. B288, 729 (1987) . 

[28] B. R. Webber, Nucl.Phys. B238 , 492 (1984) ; G. Marchesini and B. R. 
Webber, ibid. B310 , 461 (1988) . 

[29] T.D.Gottshalk , Nucl.Phys. B239, 349 (1984) ; K.Werner and P.Koch : 
CERN-TH-5575/89 . 

[30] T.Sjostrand , Comp.Phys.Comm. 27, 243 (1982) ; T.Sjostrand ,ibid. 
28 , 229 (1983) ; T.Sjostrand and M.Bengtson, ibid. 43 , 367 (1987) ; 
M.Bengtson and T.Sjostrand , Phys.Lett. 185B , 435 (1987) . 

[31] J.L.Siegrist et al. , Phys.Rev. D26 , 969 (1982) . 

[32] G.Hanson et al. ,Phys.Rev. D26 , 991 (1982) . 

[33] T. Sjostrand and Maria van Zijl , Phys.Rev. D 36, 2019 (1987) 

[34] S. Nussinov , Phys.Rev.D 50 , 3167 (1994) . 



19 



[35] V. Lugovoi , in preparation for Comp.Phys.Comm. 

[36] G.J.Alncr et al. (UA5 Collaboration) , Phys.Rep. 154 (1987) 247 

[37] C.Albajar et al. (UAl Collaboration) , Nucl.Phys. B335 , (1990) 261 . 

[38] F.Abe et al. (CDF Collaboration) , Phys.Rev.Lett. 61 , (1988) 1819 ; 
V.Barnes et al. (CDF Collaboration) , Nucl.Phys. A498 , (1989) 193 . 



20 



Figure Captions 



Fig.l In the Flux Tube model : the portion of a meson consisting of 
a segment of flux tube from the center of momentum to the i*^ quark. 

Fig. 2 The diagram of our model for the cascade break of the primary 
string A . 1 and 2 are the partons stretching the secondary string of mass M 
. Ml and M2 are the masses of daughter strings produced at the break of the 
mother string. 1 is antiquark, 2 is quark. Here and below the double line is 
diquark. 



Fig. 3 Shemes of the string breaking in our model. 

Fig. 4 a) An elementary Monte Carlo process of the Models [IS] , [2C 



: q q' + h . b) The diagram of the Lund model (and model pO|) of 
breaking of the colour singlet q'q' string into hadrons h and secondary small 
mass q"q" string , which is been brouken into two hadrons. 

Fig. 5 At the low energy Ec.m. by e^e~ annihilation a single virtual photon 
7 produces a quark- antiquark (qq) pair, tensing the primary string A , which 
decays (see Fig. 2) into secondary hadrons according to algorithm of cascade 
model of string breaking. 



Fig. 6 Illustration |^ of a hadronic event from e^e annihilation showing 



the jet axis and the components of the momentum p of a particle parallel to ( 
p\\ ) and perpendicular to ( p_L ) the jet axis. 

Fig. 7 Charged particle multiplicity distribution. — > tt+tt^ decays are 



included |^T[ , [0 . Here and below the line is the calculation according to the 
model. 



Fig. 8 Particle-density distribution a ^da/dy vs y for jets (see Fig. 6) with 
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Xmax > 0.3. Xmax IS the highest-x particle on one side of the event. The jet 
direction is oriented so that Xmax is at positive y. y is the rapidity of the 
particle relative to the jet direction assuming a pion mass. Xmax is at positive 
y and is not plotted. The distributions are normalized to the cross sections for 
jets with Xmax > 0.3. Experimental data are from pl| , p2[]. 



Fig. 9 Particle-density distribution a'^da/dpx vs pr for particles opposite 
( negative x\\ ) jets with Xmax > 0.3. pt is the component of particle momentum 
perpendicular to the jet direction (see Fig.6). Experimental data are from 



Fig. 10. The theoretical MC dependence of the average transverse mo- 
mentum of secondary strings on the mass of primary string A (see Figs. 2, 5 

)■ 



Fig.ll. The dependence of the inclusive production cross section of 
charged particles on the transverse momentum px ■ Experimental data are 
from . 



Fig. 12. The same as in Fig.ll . Experimental data are from |37 



Fig. 13. The same as in Fig.ll . Experimental data are from ||38[] . 
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